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В рамках приближения Генри (малый поверхностный потенциал и про- 
извольная толщина двойного слоя) вычислены электрофоретические подвиж- 
ности коллектива сфер для случаев постоянного заряда и постоянного по- 
тенциала. Подвижности в этих случаях существенно различны, если объем- 
ная доля дисперсии конечна и двойные слои соседних частиц сильно пере- 
крыты. Вычисления проведены на основе ячеечной модели, учитывающей 
гидродинамическое и электрическое взаимодействие сфер. 


$ 1. Введение. Коллективный электрофорез является перспективным 
методом исследования дисперсных систем [1—3]. Изучение зависимости 
скорости электрофореза от объемной доли частиц в дисперсии позволя- 
ет получить дополнительную по сравнению с электрофорезом единич- 
ной частицы информацию о двойном слое (ДС). Например, если толщи- 
на ДС х-'! велика по сравнению с радиусом частицы а 


ха< 1, (1) 

то скорость электрофореза И.; единичной сферы равна [4, 5] 
Ре 2 
Ча бл (@ 


› 


где , ц— диэлектрическая проницаемость и вязкость дисперсионной 
среды, & — потенциал поверхности частицы, Е — электрическое поле. 
И., определяется потенциалом поверхности независимо от режима фор- 
мирования ДС (постоянный заряд частицы или постоянный потенциал). 
Если же объемная доля частиц & конечна и двойные слои соседних ча- 
стиц сильно перекрыты, т. е. 


х(6—а) <! (3) 


где {=0“; 6 =а/\ — радиус ячейки, то скорость очень чувствительна, 
как будет показано ниже, к режиму формирования ДС. 

Для выявления дополнительных возможностей, заключенных в мето- 
де коллективного электрофореза, и для интерпретации эксперименталь- 
ных данных (даже качественной, не говоря уже о количественной) не- 
обходимо развивать теорию коллективного электрофореза. 

Предложенные ранее теоретические формулы {6—8] страдают опре- 
деленными недостатками. В работе Меллера, Ван Оса и Овербека [6} 
вся зависимость эффекта от объемной доли связана с изменением 
свойств равновесного ДС вследствие перекрытия ДС соседних частини. 
Этот фактор является весьма существенным, но не единственным: необ- 
ходимо учитывать еще взаимодействие электрических и гидродинамиче- 
ских полей соседних частиц. Это было сделано (как и в [6], в рамках 
модели сферических ячеек) в работах Левина [7, 8]. 
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К сожалению, в этих работах использовалось ошибочное определе- 
ние макроскопического градиента потенциала Ф, 


ЕЕ 
с0$ 9 дг 


(4) 


ГЫ 


где ф, — потенциал внутри ячейки. Ошибочность этого определения 
станет очевидной, если попытаться вычислить электропроводность дис- 
персии. Так, макроскопический электрический ток Г, равен 

К 


Де | (5) 
с05 0 07 ь 


где К — электропроводность дисперсионной среды. Электропроводность 
дисперсии Г.» равна 


ГВ 
== =К (6) 
независимо от объемной доли. Использование определения (4) привело 
к ошибке в конечной формуле для электрофореза, вследствие которой 
отсутствует предельный переход к формуле Смолуховского (по теории 
[8] скорость электрофореза сфер с тонким ДС зависит от объемной доли 
частиц). Эта ошибка компенсируется в случае использования непра- 
вильного выражения (6) для электропроводности, но на основе форму- 
лы (4) нельзя построить последовательную теорию, учитывающую, на- 
пример, поляризацию ДС. 

Правильным определением макроскопического градиента является, 
согласно [9}, 

Ф1 | -ь 


ь Ь соз 0 ® 
Его применение устраняет отмеченные противоречия. 

В данной работе использован сформулированный ранее [9, 10] на ос- 
нове термодинамики необратимых процессов метод применения ячееч- 
ной модели к анализу процессов переноса. В этих работах получены 
формулы (например, (5), (7)), связывающие макроскопические величи- 
ны с локальными полями в ячейке. | 

Использование этих определений макроскопических потоков и сил 
позволяет на основе принципа симметрии кинетических коэффициентов 
использовать формулы коллективного электрофореза для описания се- 
диментационного потенциала блока сфер. Если же применять ячеечную 
модель к рассмотрению переноса через пористые диафрагмы, то эти же 
формулы описывают электроосмос, ток течения и потенциал течения. 
Здесь мы вычислим ток течения, через который можно выразить [5, 9] 
все перечисленные выше эффекты. 

Считая б-потенциал частиц достаточно малым, мы можем прене- 
бречь влиянием возникающего при течении электрического поля на рас- 
пределение скорости. Аналогичное приближение было использовано в 
теории электрофореза Генри [4, 5]. Как и в теории Генри, будем счи- 
тать параметр ха произвольным, но не будем предполагать малости 
объемной доли частиц. 

$ 2. Решение уравнений для потоков. В соответствии с поставленной 
задачей будем считать скорость жидкости У заданной функцией коор- 
динат (она определяется из решения гидродинамической задачи без 
учета электрических объемных сил; необходимые для этого уравнения 
и краевые условия выписаны в [9], явный вид см. Приложение 1). Что- 
бы найти распределения потенциала ф, и концентрации с(г) (они опре- 
делены формулами (9,19), (9,20) !), необходимо использовать уравнения 


` Ссылки на формулы работы [19] даются в виде (9, 11) и (9, 19). 
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непрерывности для потоков ионов (9,9) 


Ри 
2= 


4 [25 (Е 2=Ф,) (- Ус (т) == Р*с(х) Уф - ра у) =0 (8) 


($ = г — безразмерные потенциалы; ШО+, 2= — коэффициенты диффу- 


зии и валентности катионов и анионов). Предположив, что относитель- 
ный перепад концентрации по ячейке мал (он связан в нашем прибли- 
жении только с несимметричностью электролита), будем искать. с(г) в 
виде ры 

с(г) ==, (1+с»(г)) (9) 


(с, — концентрация равновесного с дисперсией электролита). Локаль- 
но-равновесный потенциал ф, имеет вид 


ф = Фев + Фр (10) 


где ф., — равновесный потенциал ДС (явный вид см. Приложение 1), 
Ф» — локально-равновесный поляризационный потенциал. Величины ф,, 
бр, ф» возникают под действием внешних полей (в данном случае тече- 
ния жидкости) и линейны по ним. Если рассматривать только линейные 
явления переноса, то явный вид потенциала фр не понадобится (так как 
учет его в (8) даст квадратичные по скорости члены). 

Используя условие малости б-потенциала и формулы (9), (10), пре- 
образуем (8), выписав линейный по & член 


— — 


Д (ср + 2=Ф,) = = 1-Е (Фен) (11) 


Решения уравнений (11) имеют вид 


: 
бр 2$, = соз 6 (А* -- Ва) Е Ге | (- (=). м (дах (12) 
а 


Для определения четырех постоянных Д*, В+ воспользуемся ‚краевы- 
ми условиями (9,13), выражающими отсутствие нормальных потоков 
ионов на поверхности частицы, и двумя макроскопическими условиями. 
Согласно уравнениям переноса (9,31), ток течения надо измерять при 
отсутствии макроскопических перепадов потенциала и концентрации: 


2 


С использованием формул (9,29), (9,30) система краевых условий 
приобретает вид 


(бра ф,) |-ь=0 (13) 


д — — 
5 = 2$). =0 (14) 

Решив систему линейных уравнений (13), (14), воспользовавшись 
формулами (9,23), (9,25) для электрического тока /, и (17) (см. При- 
ложение Г) для макроскопического градиента давления Ф., получим, 
что кинетический коэффициент Г», равный 


выражается формулой 


; Ь 
и Е. 9) аз Ну | а 
те ла, 6, 9 Н (9) 4 1+2) Че (1 Р+ (5 о, | 4х ‘ ) 


(обозначения О, В, С, см. приложение [). 
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$ 3. Обсуждение результатов. Будем выражать электрофоретическую 
подвижность —Г[., в единицах =АТ/блчЕ [5, 6] и введем безразмерную 
подвижность 
6бляР 
гАТ 


0.; =_=: [4 


Интеграл в (15) можно выразить через элементарные функции и инте- 
гральную экспоненту, после чего формула приобретает вид, аналогич- 
ный формуле Генри, но намного более громоздкий. Поэтому удобнее 
предельные случаи анализировать при помощи (15), а при произволь- 
ных {| и ха пользоваться численным методом. 

Главным (но не единственным) фактором, определяющим зависи- 
мость П., от объемной доли, является перекрытие ДС соседних частиц. 
Чтобы сделать анализ более наглядным, разделим всю область измене- 
ния параметров ха и \ на две части (рис. 1). Будем условно считать, что 
при х(Б—а)< 1 (область 1) ДС существенно перекрыты, а при 
х(5—а)>1 (область 2) —слабо перекрыты. В области 2 (по крайней 
мере, достаточно далеко от границы) справедлива формула Генри и 
учет взаимодействия частиц дает лишь небольшую поправку (последнее 
связано с тем, что в рассматриваемом приближении поляризация ДС не 
учитывается. С учетом поляризации ДС взаимодействие частиц может 
дать значительный вклад в эффект и при ха»1). Погрешность форму- 
лы Генри по сравнению с более общей формулой (15) растет по мере 
приближения к границе областей и становится очень большой в обла- 
сти 1 (рис. 2). При ха=. (предельно толстый ДС) 0.,=0 при постоян- 
ном б-потенциале и равна 
й. 3. 9 (5. 

. 2 НС — 1) 


при постоянном заряде (обозначение &.. см. приложение 11) *. Такое 
поведение 0., можно пояснить при помощи наглядных соображений. 

При *б—0 потенциал ф.г и объемная плотность заряда в ДС р не за- 
висят от координат. Поэтому ток течения равен 


(16) 


Ь—о[ 
где /, — объемная скорость жидкости, а 


= — Роб (+ г) при 2: сопзЕ 


© = — т при ву == сопз${ 
Ап (52 — аз) 62 (1— 7?) у 


3 


При Е==соп${ и хб->0 (что соответствует с.>0) плотность поверхност- 
ного заряда частицы и плотность объемного заряда стремятся к нулю. 
При *Ь==0 частицы вовсе незаряжены и поэтому не движутся в электри- 
ческом поле (0.,==0). При постоянном д, плотность объемного заряда 
не зависит от х, но убывает с ростом радиуса ячейки как 65-*. Однако 
скорость /, при заданном внешнем давлении растет пропорционально 
5, поэтому ток течения в этом предельном случае конечен. 

Как видим, при сильном перекрытии ДС электрофорез становится 
весьма чувствительным к режиму формирования ДС, поэтому коллек- 
тивный электрофорез может стать дополнительным источником инфор- 
мации`о ДС. 


? Выражения «постоянный б-потенциал» и «постоянный заряд» означают, что эти 
величины не меняются с изменением у и ха в любых пределах. 
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Рис. 1 


Рис. 1. Области параметров, при которых сущест- 
венно (1) и несущественно (2) перекрытие ДС 
соседних частиц 


Рис. 2. Отклонения от формулы Генри при конеч- 

ной объемной доле; ха=1; 1 — формула Генри; 

2 — формула (15) при 060—015 8 — формула 
(15) при &=с01п8+ 


Рис. 3. Влияние гидродинамического взаимодейст- 

вия частиц на электрофорез; б==соп$ё; 1— 

ха-=0, модель Хаппеля; 2 — ха=0, модель Ку- 

вабары; 3 — ха=0, модель Мехта и Морзе; 4 — 
ха=1, модель Кувабары 


Рис. 3 


На рис. 3 приведены графики зависимости (16) для различных гид- 
родинамических ячеечных моделей. Различие между ними не превыша- 
ет 15% при всех у. Интересно отметить, что множитель Н(\), введен- 
ный [9] для других целей (выполнения формулы Смолуховского), улуч- 
шает согласие результатов различных ячеечных моделей (без него раз- 
личие превышало 50%). | 

Предельный случай (16). интересен еще и тем, что показывает влия- 
ние взаимодействия частиц на электрофорез «в чистом виде» (толщина 
ДС и степень перекрытия не влияют на эффект). Все убывание й., с ро- 
стом 4 (кривая 2, рис. 3) обусловлено гидродинамическим взаимодей- 
ствием и не может быть получено в рамках теории [6], не учитывающей 
этого взаимодействия. Поэтому вывод. авторов работы {10 том, что 
наблюдавшаяся в их опытах независимость й.; от у служит подтвержде- 
нием теории [6], представляется необоснованным. Кроме того, условия 
эксперимента в [1] не дают возможности количественной интерпретации. 
данных на основе Теории [6] или предлагаемой теории из-за высокого 
5-потенциала (89 мВ) и полидисперсности системы. 


Приложение 1. 


Явный вид поля скорости. Поле скорости в ячеечной модели имеет 


вид [11]: 
о, = —- [с 39. (5) +2 (- } а, (+) в 
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{нам понадобится только радиальная компонента). 
Макроскопическая объемная скорость равна [9] 


= 0/Н (1) 
Макроскопический градиент давления равен [9] 


Ф,= За Н (у) С. 
в 9 (1) 


(17) 


Величины @, Н, С, —С. для различных ячеечных моделей 


ОЗ ЕЕ ЕЕ ОБЬ ЕВА 


Величина Модель Кувабары Модель Хаппеля Модель Мехта и Морзе 
А ООО ООО 
С, 1 18/2 1--3у5/2 "1 - 53/4 — 55/4 
С» —1/4 — 38/40 —1/4 —(1 — 1 55) /4 
Сз я —5 —15( — 19/2 
Са —3/2 —3(4--25/3) /2 —3/2 
9(7) 1—9%/5 -- уз — 16/5 1—3%/2 + 35/2 — 8 |1 — 9/4 + 53/2 — 55/4 — 
— 8/2 8/2 
1—5 1—53/2 -|- 55/2 
НО 1 Е. дат 94 
— Ч—9 +=) у 


Приложение И. 


Явный вид равновесного потенциала ДС. Приводимое здесь решение 
более подробно обсуждается в [8]. Уравнение Пуассона — Больцмана 
при малых &-потенциалах имеет вид 


га а- — 
а” Ча (18) 


‘где х=хг; х==ха, х.=иб; краевые условия имеют вид 


Фе „= 5 
а- 
ях РеЕ ные . 


(последнее краевое условие обеспечивает электронейтральность ячейки 
как целого). 
Решение (18) имеет вид 


Фев — С (Ра ехр (х) + Рьехр (-— 4))/х 
Р.= (х, +1) ехр (—х,)/Р; 
Р»= (х,—1) ехр (х.) /Р, 
Р;= (+1) ехр (ж—х!) + (х.—1) ехр (х.—х) 
Плотность поверхностного заряда в, связана с & формулами 


Е — _4подраР 
ы гЮТ 


бе — Л -- Хо) 
= ыЫРз/(ж - 1) (<, — Юехр(х, — ж) — щ%—1) (-Е1) ехр (ж, — х,)) (19) 


5 имеет смысл потенциала единичной сферы без ДС (х,=0); &, — по- 


тенциал единичной сферы при х50; 5 — потенциал сферы в ячейке. 
Для получения предела (16) удобно использовать предельное значение 


{18) 
Зв? 


"ва-м при ж—0; у =-0 
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агЬНгагу 1сКпез$ оЁ {Пе доцЫе 1ауег) Пе е!есгорНогейс то ИНез оЁ Ше соЦееНуе о 
зрНегез Вауе Бееп са1сша{е@ {ог \1е сазез оЁ а сопзфап{ сВагое ап@ а сопфапё роепНа]. 
[{ Баз Бееп зНо\ми 4На{ 1 \10о$е сазез {!е тоИЙез аге зиб$фапНаПу-. Чегет, И Ше 
ЪиК рогНоп о{ @1зрегзюп 1$ НийЙе ап@ Ше доцЫе 1ауегз о! а@]асеп{ рагИс!ез аге тагКе4- 
1у оуеПаррей. Са!сфаНопз Бауе Бееп сагуе@ ошё оп Фе Баз1$ оЁ а сейаг тоде! {аки 
Ио ассоит{ {Не Будтодупап с ап {Пе @есёса! ищегасНоп оЁ зрНегез. 


